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ABSTRACT 

The gamma-ray burst (GRB) jet powers the afterglow emission by shocking the surTounding 
medium, and radio afterglow can now be routinely observed to almost a year after the explo¬ 
sion. Long-duration GRBs are accompanied by supernovae (SNe) that typically contain much 
more energy than the GRB jet. Here we consider the fact that the SN blast wave will also 
produce its own afterglow (supernova remnant emission), which will peak at much later time 
(since it is non-relativistic), when the SN blast wave transitions from a coasting phase to a 
decelerating Sedov-Taylor phase. We predict that this component will peak generally a few 
tens of years after the explosion and it will outshine the GRB powered afterglow well-before 
its peak emission. In the case of GRB 030329, where the external density is constrained by 
the ^ 10-year coverage of the radio GRB afterglow, the radio emission is predicted to start 
rising over the next decade and to continue to increase for the following decades up to a level 
of ^ mJy. Detection of the SN-powered radio emission will greatly advance our knowledge 
of particle acceleration in ^ 0.1c shocks. 

Key words: radiation mechanisms: non-thermal - methods: analytical - gamma-ray burst: 
general - supernovae: general 


1 INTRODUCTION 


The external shock, produced as the relativistic gamma-ray hurst 
(GRB) jet interacts with the circumhurst me dium, is thought to give 
rise to the GRB afterglow radiation (e.g., |Sari. Pir an & Narav^ 
1 19981 : IWiiers & Galarndll999l : IPanaitescu & Kumaj 20QO|) . Radio 
afterglows have been observed for up to a decade after the time of 
the explosion, when the GRB jet blast wave moves at subrelativistic 
speed. The total kinetic energy in the blast wave i n GRBs is 10^ ^ 
erg (after beaming correction is applied, see, e.g.. lFrail et al.l200lh . 

Long-duration GRBs are associated with supemovae (SNe) 
of the rare broad-line Ic type (see, e.g, Wooslev & Bloon] 


__ ;.g. IWoc 

l2006l : !Hiorth & Bloomll2012l : lMelandri et al. Il20l4 and references 

therein). These peculiar SNe have very large kinetic energies, 
which are 50 times larger than the total (beaming-corrected) 
GRB energy, and their ejecta have also unusually high velocities 
(almost ~ 0.1c). The SN ejecta will also interact with the external 
medium: it will also drive an external shock, which will produce 
an afterglow. This SN afterglow is simply the emission from the 
supernova remnant (SNR), therefore, we will refer to it as “SNR 
emission.” Since the total energy in the SN is much larger than the 
one in the GRB jet, the SNR emission will eventually outshine that 
of the GRB. 

In this work, we calculate the expected radio synchrotron 
emission from the SN external shock. Due to the non-relativistic na- 
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ture of the SN ejecta, it will only decelerate after a few tens of years. 
As in the case of the GRB jet, the synchrotron emission from the 
SN blast wave depends on the external density and the shock micro¬ 
physics, which are unknown quantities a priori', however, we will 
use our knowledge of GRB afterglow studies as a guideline. In par¬ 
ticular, GRB 030329 has been followed up in the radio for a decade 
after the burst and its light curve continues on a smooth, power-law 
decline implying constant parti cle acceleration efficiency when the 
shock becomes subrelativistic jMesler & Pihlstr(inill2013l) . We use 
the available sample of GRBs associated with SNe to predict the 
SNR radio emission. We find that, within our assumptions, the SNR 
radio emission should be detectable, and lead to a rise of the radio 
light curves about a decade after the GRB trigger. 


There are two long GRBs at low-enough redshift to allow for 
the optical detection of the accompanying SN at the same level 
of SNe accompanying GRBs (GRB-SNe). However, for these two 
bursts no SN optical emission was detected l Della Valle et alj2006l: 


iFvnbo et al.l2()06l : lGal-Yam et al.l2Q06l : lGehrels et alj2006h . These 

two cases imply that either no SN explosion took place or that the 
SN was optically faint, because, e.g., of a very low amount of ®®Ni 
present in these explosions. If the kinetic energy of the SNe of these 
two bursts was similar to that of other GRB-SNe, then we can use 
our model to predict their radio SNR emission. This radio emis¬ 
sion, if detected, would imply that a SN of the GRB-SNe type was 
present in these two GRBs. 

This paper is organized as follows. In Section we summa¬ 
rize the model we adopt for the synchrotron emission from the ex- 
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temal shock of the GRB jet in the Sedov-Taylor regime. In the same 
Section, we calculate the synchrotron emission expected from the 
external shock of the SN ejecta. In Section]^ we take the known 
sample of GRBs with a SN association and calculate the expected 
SNR radio emission. We do the same with the two GRBs with no 
SN association mentioned above. We finish with a discussion in 
Section|4] 


2 THE RADIO AFTERGLOW AT LATE TIMES 


The GRB jet interacts with the external medium and drives a 
beamed blast wave that accumulates mass and decelerates. As it 
reaches non-relativistic velocities, the b last wave ent ers the Sedov- 
von Neumann-Taylor (ST) phase (e.g.. lTavloilll946ll and becomes 
approximately spherical (albeit with a displaced center from the 
center of the explosion, see Section [13. Since we are interested 
in late time radio observations, we focus our discussion on the 
synchrotron emission i n this phase. We present a short summary 
following the wo rk of Sironi__&_Giannioj (|2013|) (see references 
therein, and. e.g..|Chevalieil 1982al f5. Il998t Li & Chevalieilll99^ : 


IChevalier & Franssonl200d . for a similar modeling of the dynamics 

and radio synchrotron emission of the SN blast wave as it interacts 
with the surrounding medium). 


2.1 The GRB jet afterglow 


For the GRB jet compon e nt, a cc ording to numerical 
studies (e.g., [Wygoda et al.| 1201 ll : IPe Colle et al.l |2012| : 
Ivan Eerten & MacFadvenI 20121) . the spherical ST solution 
can be used to estimate the observed flux after a time tsx, given by 


tsT ~ 290(i5GRB,5l/tlo)^^^(l + z) d. 


( 1 ) 


Beyond tsT, the blast wave radius, R, and velocity, /3 (in units of 
the speed of light), are 

R « 9.4 X 10^^(EGRB,5i/no)^^®fy^®(l + z)“^^®cm, (2) 

P « 0.4(i5GRB,5i/no)^^^ty/'^®(l + (3) 


where Egrb is the beaming-corrected kinetic energy of the GRB 
jet blast wave, n is the number density of the external medium, 
t is the observed time since the explosion, z is the cosmolog¬ 
ical redshift, and we have used the common notation Qx = 
<3/10“^ in c.g.s units. We assume here a constant density external 
medium, although a wind medium is also a possibility (see, e.g., 
IChevalier & Li|[l99^ . At this stage, the synchrotron emission for 
an observed frequency v (where < u < Uc, and Ua, 

Vm and Vc are the synchrotron self-absorption, minimum injection 
and cooling frequencies, respectively) is 


F. 


X 


(0.02 mJyje^JL/ 


1+p 3+5p 19-5p 

^ 4 p 10 n 20 

^S,-2-^GRB,51^^0 


3(7-5p) 1-p 5p-8 o 


(4) 


where Ce = 4ee(p — 2)/(p — 1), te (es) is the fraction of shocked 
fluid energy in electrons (magnetic field), p is the power-law index 
of the energy distribution of the electrons and d is the luminosity 
distance. The normalization corresponds to p = 2.4, but it does not 
depend strongly on p. We assume and es do not change with 

time. _ 

As discussed in ISironi & Giannio j j2013l) . caution must be 
taken when the minimum Lorentz factor of the electrons, 'ymin, 
given by 'ymin — 1 = (l/8)ee(m-p/me)/l^, drops below ~ 2. At 


this point, the system will transition to a new regime: the “deep 
Newtonian” (DN) regime, where the spectrum of accelerated elec¬ 
trons can no longer be approximated as power-law of index p in 
kinetic energy, but it still follows a power-law distribution in mo¬ 
mentum with slope p. In this case, for 2 < p < 3, the bulk of the 
electron energy is contrib uted by mildly rela tivistic particles with 
Lorentz factor of ~ 2 (see lGranot et al.|[2006h . This happens when 
the velocity drops below 

PuN~0.2e-li\ (5) 

which occurs at a time fDN, see eq. 10, 

tDN « 2.2{Egw,,&i/ noY^^e^L\(l +z)yr, (6) 


and the synchrotron flux in this regime is given by 


Fn 


1+P 11 + P 3+3p 

« (0.01 niJy)ee,-l£B(‘_2^'GRB,5l’^0 

-3(l+p) 1-p 4-p « 

X V ^GHz(l + ^) ^ ^27.5- 


(7) 


2.2 The SNR emission 

Long GRBs are (usually) accompanied by a powerful (and fast ) 
SN, see, e.g.. IWooslev & BloomI (l200fil) : iHiorth & BloornI ll2012h . 
As the SN propagates through the external density, a SN blast wave 
will develop, which will also produce its own afterglow. We refer 
to this signal the “SNR emission”. In essence, it is similar to the 
GRB jet afterglow, but the SN blast wave is much slower (while 
the GRB jet is relativistic), and the SN energy (at least for SNe 
which accompany GRBs) is much larger than the typical (beaming- 
corrected) energy of GRB jets (see Table[Tll. 

The geometry we have in mind is the following (see fig. 7 
of iRamirez-Ruiz & MacFad^^l2010l) . Both the GRB blast wave 
(at late times) and the SN blast wave (always) are approximated as 
spherical blasts. However, the GRB blast wave, because of its initial 
directionality, will be expanding as a (quasi-) sphere with a shifted 
center. This new center is located at a distance (from the central 
engine) at which the initially beamed GRB jet reaches a blast wave 
Lorentz factor of 2. The SN blast wave, on the other hand, is 
a truly centered (on the central engine) explosion. It probes mainly 
the pristine “unshocked” external density along the equatorial plane 
(not affected by the GRB jet blast wave) rather than the one in the 
poles. Once the SN blast wave reaches a size similar to that of the 
GRB blast, it will outshine the latter. 

The SN blast wave coasts with constant velocity until it starts 
to slow down after it has doubled its initial mass. The deceleration 
time is 

fdec,SN « 29/3gN(^®j^(iJsN,52.5/wo)^^^(l + z) yr, (8) 

where Psn and Esn are the velocity and energy of the SN blast 
wave, respectively. At this time, the SN blast wave velocity will 
decrease as oc see eq. 0. The SN blast wave radius will 

increase linearly with time while the velocity is constant, but will 
increase slowly with time, as oc after fdec.SN, see eq. 0. 

Fig.E shows the velocity of both components, the GRB and SN 
components, as a function of observer time and radius. 

During the coasting phase, while the blast wave collects more 
and more external medium, the SNR emission light curve will 
rise rapidly as oc (see below). Because of its subrelativistic 
speed /3 sn <0.1, the SNR emission is in the DN regime, see eq. 
0. For this reason, after fdec.SN, the SNR emission will decay 
as oc f“3(i+p)/io^ described in eq. 0, and the peak flux (at 


© 2015 RAS, MNRAS 000, 000-000 















































GRB radio afterglow rebrightening 3 


GRB 

SN 

Z 

V [10® km/s] 

Esn [10®! erg] 

980425 

1998bw 

0.0083 

24 

50 ±5 

030329 

2003dh 

0.1685 

29 

40± 10 

031203 

20031W 

0.106 

21 

60± 10 

060218 

2006aj 

0.0335 

19 

2 ± 0.5 

100316D 

2010bh 

0.0593 


~10 

101219B 

2010ma 

0.55 



120422A 

2012bz 

0.283 

20.5 

41 

130702A 

2013dx 

0.145 

24 

35± 10 

130427A 

2013cq 

0.34 

32 

64 ±7 

140606B 

iPTF14bfu 

0.384 




Table 1. GRB, its associated SN, redshift, the SN velocity at 10 d 
after the explosion and the SN energy. GRB-associated SNe are typi¬ 
cally very energetic, with several xlO^^ erg, and fairly fast velocity of 
~ 0.1c. In our modeling, for blan k cells in this tab l e we a ssume 1998bw- 
like values. Data from table 9. 2 in|Hjo rth & Bloonj j20 12[) (and reference s 
therein); additio n al data from 

IXu al] j2013h: iLeyan e t alj j2Q14h : ISchulze et alj i2Ql4) ; iD’Elia et alj 
i2015l) : ISinger et al.l i2015t) . 



actually provide an independent constrain on £^grb, An and den¬ 
sity (which does not depend on microphysics). This constrain de¬ 
pends most sensitively on the SN velocity. On the other hand, if we 
have a good knowledge of Eqrb and /Ssn we can determine the 
external density. We note that for t > fsN,dec, the SNR emission 
exceeds that of the GRB by a factor ^ 1, 

see eq. Q. For typical values, the flux ratio is > 100. 

The expressions presented in this Section (for both the GRB 
jet and SN components) are applicable when the radio observing 
frequency is max{Ua,i'm) < y < Vc (the cooling frequency is ir¬ 
relevant in the radio band). For the SN component, since the peak 
of the emission is produced when the blast wave is far away from 
the center of the explosion, and moves at a low velocity, it can be 
shown that bo th Vg and Um are always below the observin g ra¬ 
dio frequency jNakar & Piranll20lll : lsironi & Gianniosll2013h . The 
same is usually true for the GRB jet component emission at very 
late times. 

We note that the SN blast wave radius, i?sN, at f > fdec.SN 
increases as ^2/5 

, see eq. The angular size of the SNR emission 
at this stage would be 0 sn ~ 2Rs^/dA, where dA = d/{l + z)^ 
is the angular diameter distance. The SN blast wave radius at this 
stage only depends on iSsN, density and time, see eq. 13. We can 
solve for density and substitute it in the SN blast wave synchrotron 
flux equation [eq. 0, but with iSsN instead of Sgrb]. We find 
an interesting expression which does not explicitly depend on time 
since the explosion (although the angular size and flux do vary with 
time), as follows 


i+p 

— 4 

^e,-l ^b,-2 


30 


e 


X y. 


1 mjy / \ 1 mas 

3-1 , , ll+3p 

■2-(l I ,'|-(2+p) J 4 

:Hz U + “27.5 ■ 




( 12 ) 


Figure 1. Velocity (in units of the speed of light) as a function of ob¬ 
server time and blast wave radius. The solid (dot dashed) line marks the 
GRB jet (SN) component. The SN component has a constant velocity, until 
fdec SN. when it decelerates. These lines coiTespond to the parameters of 
GRB 030329, described in the caption of Fig.|2] At t > 20 yrs, the SN blast 
overtakes that of the GRB jet and its emission dominates. 


Strictly speaking the normalization is only valid for p « 2.4. 
Therefore, by measuring at a given time the flux and size after the 
peak of the SN radio light curve, we can constrain ~ CeCs with our 
knowledge of i?sN. 


fdec.SN) will be given by 

1 + P 1 + P 

Fjy(t = tdec,SN) ~ — _1 

^+P 1-P 1-p _2 

X F;sN, 52.6 tlQ ■* i^GHz(l + ^) ^ '^27.5- 
For t < isN,dec, the SNR radio flux is given by 

1+P ii+P 

^deCjSN) ~ 

3+P „ 1-P -(5 + p) _2 

X nQ tlOyi^GHz E z) 2 £(27.5- 


(9) 


( 10 ) 


Old, nearby and fast moving SN ejecta, which interact with a dense 
external medium, would yield larger fluxes in this stage. 

We can find a time, teq, where the declining GRB afterglow 
and the rising SNR emission flux are equal. At these late times, the 
two blasts propagate at a similar speed and at ~ pc scales away 
from the explosion center. Assuming constant density surrounding 
gas and the same microphysical parameters (ce, es, and also p) we 
find 


feq = tdec,SN{EaRB/EsiSlf^^ 

« 9,0^N.'^®i('E’GRB.5i/no)^^^(l + 2:)yr. (11) 

In view of this last equation, determining observationally feq can 


3 APPLICATION TO GRBS 
3.1 Sample of GRBs 

We present the current sample of GRBs for which there is a strong 
association with SNe, see Table[T] We provide the approximate ve¬ 
locity of the SN (at ~ 10 d) and the estimated total energy Esn- 
With this information, we can predict the SNR light curves for a 
given external density and microphysical parameters. 

Alternatively, we can model the very late radio GRB afterglow 
data (~ years time-scale) in the context of synchrotron emission 
from the GRB jet blast wave and determine the density and mi¬ 
crophysical parameters. The medium that is probed by the blast 
wave at such long time-scales (large distances), and the micro¬ 
physics of the blast wave when the velocity is small, would be a 
good proxy of the medium/microphysics also expected for the SNR 
emission. This exercise is only possible in the case of GRB 030329 
for which there is radio af terglow data for ~ 8 yrs after the burst 
dMesler & Pihlstr(imll2013t) . We will treat this case separately. In 
general, GRB afterglows are detectable for less than a year. In these 
cases, the external density probed might be very different than the 
one probed by the SN blast wave at very large distances (the same 
applies to the microphysics), and any attempts made to connect the 
early GRB radio emission with the late SNR one would be plagued 
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Figure 2. Top panel: We fit the radio 4.9 GH z light curve (black circles, 
iMesler et alJIlOlj ; iMesler & Pihlstron]|2013l and references therein) with 
the lowest possible density allowed by the data of n = 0.2 cm~^ (dashed 
blue line), see Section [52] With these parameters, we predict that the SNR 
radio emission (dot dashed blue line) should rise above the EVLA flux limit 
of ~ 20/iJy (horizontal black dashed line) by 2030 at the latest (blue solid 
thick line marks the sum of both the GRB jet and SNR emission compo¬ 
nents). The year 2015 is marked with a vertical green dashed line. Strictly, 
speaking, our fit to the radio data with eq. (4) is valid only for t > fsT 2 
yrs (same applies to the bottom panel). We only include radio data with 
t > 100 d for both panels. Bottom panel: Measured angular sizes (black 
circles) and upper li mit (triangle) of the GRB 030329 radio afterglow source 
iMesler et alfcoij and references therein). The blue dashed (dot dashed) 
thin line is the GRB (SN) afterglow angular size for the model presented 
in the top panel. The expected angular size at a given time is the maximum 
between both components. The dashed magenta vertical line marks when 
the “Deep Newtonian” regime starts for the GRB component (the SN com¬ 
ponent is always in this regime). 


by this uncertainty. For this reason, for these GRBs we will leave 
the density/microphysics as free parameters when we predict the 
SNR emission. We will also treat GRB 980425 separately, since it 
is the closest GRB observed to date. We will then consider the rest 
of GRBs in Table^for our analysis. 

In addition to the GRBs in Tabled we will also consider two 
GRBs for which a SN was expected but not found, GRB 060505 
and GRB 060614, at 2 = 0.0894 and 2 = 0.125, respectively 
jPella Valle et alJ I 2 OO 6 I: Ipvnbo et alj bood: lOal- Yam et all bood: 
iGehrels et al.ll2006l) . For these two, we will use 1998bw as a tem¬ 
plate to predict the expected SNR emission, which could serve as 
a test whether or not a SN like 1998bw was truly present in these 
two explosions. 

We note that throughout this paper we set an approximate limit 
of detection of a ra dio signal of ~ 20 /jJv, which is expected hy 
EVLA at 4.9 GHz jPerlev et alJbOl ih . We will use this limit to 
determine if and when the SNR emission will he ohservahle. 


electrons p « 2.2). We can use the same density/microphysical 
parameters used to model the radio GRB afterglow to predict 
the SNR emission, since at such long observed time-scales both 
the GRB and the SN components move at similar velocities and 
probe the density of the external medium at similar scales. Sev¬ 
eral groups have fitted the radio afterglow (see, e.g.. IMesler et al.l 
b012l : lMesler & Pihlstroml2013l . and references therein), and a few 
sets of density/microphysical parameters have been provided. Fits 
point out to a uniform external medium of ~ 1 cm“®. 

As can be seen from eqs. l[^ and l|3l, the larger the density, the 
earlier and brighter the SNR emission peak will be. For this GRB 
we take a conservative approach and estimate the lowest density 
allowed by the data. We do this hy setting the maximum values 
of the microphysical parameters (ce = es = 1/3) and we set the 
energy to the larg est value found for this GRB of 1.5 x 10®^ erg, see 
tables 2 and 3 of IMesler & Pihlstrdml (l2013l) . For a fixed observed 
flux, these choices will yield the minimum required density, see eq. 
(0, which is n « 0.2 cm This constrains the rise of the SNR 
radio emission to occur on or before ~2030, see Fig. We also 
include the afterglow size (both the GRB and SN components) for 
the same model. As can be seen, the SNR emission size dominates 
the GRB afterglow size ~ 20 yr after the explosion. 

Using the last observed radio point at ~ 8 yr as a limit of when 
the rise of the SNR emission occurred, e.g., feq > 8 yr, we can find 
an upper limit on the external density (for a given /3 sn and Egrb), 
see eq. GD- For the parameters mentioned above, we find that the 
density cannot be larger than 4 cm~®, otherwise we would have 
seen the rise of the SNR emission already. 


3.3 GRB 980425 


We now focus on GRB 980425. For this GRB, there is no radio 
data after about a year, therefore, it is not possible to follow the 
same procedure as for GRB 030329. To calculate the SNR radio 
emission for this burst we simply assume a typical value found for 
GRB afterglows Ce ~ 0.2, and an optimistic value of es ~ 0-01 
teamiol Duranlb014 : Santana. Barniol Dura n & KumailbOldl see 
also simulations by Sironi & St)itkovskv(l201 1 ). We also use p — 
2.4. We use these values of Ce, es and p throughout the rest of the 
paper, unless otherwise noted. We allow for two possible values for 
the external density, n = 0.1 cm~® and a more optimistic value of 
n = 1 cm“®. Using the velocity and energy of the SN component, 
we can predict the SNR emission, and it is shown in Fig. |3 As 
mentioned above, a higher (lower) density yields a SNR emission 
that peaks earlier (later) and is brighter (weaker), but the shape of 
the afterglow is preserved. For the parameters adopted for this SNR 
emission, it will reach values at the same level of the GRB radio 
afterglow detected < 1 yr after the explosion. 

We use GRB 980425 as a benchmark case. Assuming an opti¬ 
mistic density of 1 cm“®, we find that a GRB with a 1998bw-like 
SN should be located at a distance < 1 Gpc (z < 0.25), so that its 
SNR radio emission peaks above ~ 20/rJy. 


3.2 GRB 030329 

GRB 030329 holds the record for the longest radio afterglow 
ever detected. The last reported data point is at ~ 8 yr 

dMesler & Pihlstr(imll2013l) . We can fit the late radio data, which 
decays as , using the synchrotron flux of the GRB jet blast 

wave, see eq. 0 (which points out to a power-law index of 


3.4 Rest of the sample 

For the rest of the GRBs in Table [T] we follow a similar proce¬ 
dure as the one followed for GRB980425. For GRB 101219B and 
140606B, we assume a 1998bw-like SN. For GRB 100316D, we 
assume a 1998w-like SN velocity. We predict their SNR radio emis¬ 
sion and show them in Fig. |4l along with their available late time 
radio data. We also include the two GRBs with no SN identification 
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Figure 3. it Top panel: The 4.8 GHz radio observations of GRB 98 0425 are 
denoted with black circles jKulkami et al.|[l9^ : lFrail et al.ll2QQ3h . We use 
eg = 0.2, es = 0.01 andp = 2.4 to calculate the SNR emission, with n = 
1 cm~^ and n = 0.1 cm“^ (dashed and solid blue lines, respectively). The 
vertical green dashed line marks the year 2015. Bottom panel: The size of 
the SN blast wave for the same densities presented above; solid and dashed 
lines overlap before 50 yrs, since the coasting velocity is kept fixed, and 
only the density is changed. 


(GRB 060505 and 060614) assuming that a SN blast wave with the 
properties of 1998bw was present. For the parameters adopted in 
Fig. a it can be seen that for some bursts the SNR radio emission 
could potentially be detected, while for others it is too weak to be 
detected with current instruments. 

As the external density is an uncertain parameter, we allow 
for a larger (smaller) density than the adopted value of 1 cm~® by 
a factor of 10. We indicate the position of the peak SNR emission 
for these different densities with a dashed arrow in Fig. 0 We do 
this only for the case of GRB 100316D, which gives a sense of the 
density dependence of our calculation. 


3.5 Host galaxy radio contribution 

The host galaxy contribution to the observing radio band is a criti¬ 
cal factor to be able to identify the SNR radio emission a few tens 
of years after the burst. For example, the host galaxy contribution 
of GRB 980 425 appears to be bright, on the order of 420 ± 50 /rJy 
at 4.8 GHz dMichalowski et al.ll2009l) . Also, the host galaxy con¬ 
tribution of GRB 031203 is estimated to be 216 ± 50 /tJy at 5.5 
GHz dStanwav, Davies & LevanI 1201 OL see, also, ISoderberg et al.l 
l2004h . Other bursts in our sample have a host contribution in ra¬ 
dio which is smaller: 3 -(t limits at 5.5 GHz of 117, 37 and 33 
^Jy for GRB 060218, GRB 060505, GRB 0606014, resp ectively 
dStanwav. Davies & Levanifioid i lMichalowski et al.ll2012l) . There¬ 
fore, setting a precise baseline observation to reduce the errors in 
the host contribution measurements is of utmost importance to be 
able to detect the radio SNR emission. 



Figure 4. The left-hand panel shows the observed radio data (circles) and 
upper limits (triangles) as a function of time since the trigger (in log- 
scale) for some bursts in our sample (see legend for colors). The right- 
hand panel shows the predicted SNR emission as a function of absolute 
observing time in years (in linear scale). For these light curves we use 
Ce = 0.2, tB = 0.01, p = 2.4, and an external constant density of 
n = 1 cm“®. The velocity and energy of the SN component can be 
found in Table [T] and are used to obtain the SNR emission light curves. 
For a higher (lower) density of n = 10 cm~® (n = 0.1 cm~®), the 
SNR emission will peak at an earlier and higher (later and weaker) flux, 
indicated by the blue dashed arrow (see GRB 100316D as an example). 
The horizontal black dashed line marks the EVLA flux limit of ~ 20;rJy. 
For GRB 060505, and 060614, we assume a 1998bw-like SN blast wave. 
Radio data (prese nted only for some bur sts) and predicti ons at: 4.9 GHz 
for GRB 032103 {Soderber^ et all l2004h. GRB 0 60218 JSoderberg et alJ 
l2006j). GRB 060614 jLondish. Wieringa & Fraill l2006h GRB 130427A 
Ivan der Horst et al.ll2014h. GRB 1.30702A Ivan der Horslll2013ll : 5.4 GHz 
for GRB 1003160 jMarguttietal.l 12011 : 6.1 GHz for GRB 14 0606B 
ISinger et a i]l201.5h: 8.5 GHz for GRB 060505 iOfek et al.ll2007h GRB 
101219B <Fraii201ll) . and 15 GHz for GRB 120422A jSchulze et alj20T4l) . 


4 DISCUSSION 

As the GRB jet interacts with the external medium, it produces an 
external shock, which gives rise to the observed GRB afterglow 
emission. The SN ejecta, which accompanies long GRBs, will also 
interact with the external medium and will also drive an external 
shock. We have calculated the expected synchrotron emission from 
this shock and predicted the early phases of the radio SNR emission 
for the current sample of GRB associated with SNe. 

The prospects for detection depend on the assumed values of 
the external medium and microphysics of the SN blast wave, which 
are a priori unknown and are the free parameters in our model. At 
least for the case of GRB 030329, where the GRB radio afterglow 
was detected to about ~ 10 yrs after the explosion, we have esti¬ 
mates for the external density at large distances from the explosion 
(parsec scale), and also on the microphysics of a blast wave that 
moves close to ~ 0.1c, and we can use these to calculate the SNR 
emission. For this burst, the prospects for detection are optimistic. 
However, if the SN ejecta shock microphysical values are much dif¬ 
ferent than those for the GRB shock, then the expected SNR radio 
emission would be much harder to predict, since the peak of the 
SNR emission light curve depends linearly (almost linearly) on the 
fraction of shocked energy that goes into electrons (magnetic field). 

The value of the power-law index of electron energy dis- 
tribution expected in non -relativistic shocks is p « 2 (e.g., 
iBlandford & Eichl^l 19871) : however, we use a slightly larger con¬ 
servative value of p « 2.2 — 2.4. We also assume that Ce is of 
the order of 10 per cent, guided by GRB afterglow studies and by 
the GRB afterglow fits done for GRB 030329. It is possible that 
Ee is smaller for a non-relativistic blast wave, which will weaken 
the SNR radio signal significantly (linearly, as mentioned before). 
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Since the magnitude o f £e for subrelativistic shocks contin ues to he 
actively studied (e.e. JPark. Canrioli & SDitkovskvll2014h . we take 
the value adopted in this paper as an optimistic one. Nevertheless, 
SNe are well-known active sites of particle acceleration and radio 
emission. In fact, mildly relativistic SNe, which radiate in the ra¬ 
dio band, have been detected even without the presence of a GRB; 
“typical” Ibc SNe efficiently acceler ate particles in ^ 0.1c shocks, 
which produce radio emission ('e.g.. lSoderberg et alj|20ld) . Thus, 
expecting a radio signal from SNe accompanying GRBs is a natural 
extension of what we presently know, and would provide an excel¬ 
lent opportunity to study the behavior of Ce as the shock transitions 
to ~ 0.1c. 

The radio SN rebrightening and the large size of the SN blast 
wave provide a unique opportunity to directly resolve these sources, 
if the GRB is located al z< 0.2. In the coming years, GRB 030329 
and 980425 may reach flux levels in excess of 100 /rjy and source 
size of > 1 mas 20 mas for GRB 980425 / SN 1998bw), well- 
within the resolving capabilities of current instruments. 

We have also assumed that the external density probed by the 
SN blast wave, at ~ several pc from the center of the explosion, 
is constant and roughly ~ 1 cm~® (although we have considered 
also values larger and smaller than this by a factor of 10). To ex¬ 
plore these assumptions, let us assume a GRB progenitor Wolf- 
Rayet star with wind velocity of ~ 1000 km/s and mass loss rate 
of ~ 1O“®M0 yr“^ that extends to infinity. At parsec scales, the 
wind density drops to ~ 0.03 cm~®. It is safe to assume that at 
some point, this density will encounter a “floor.” If the star is in 
isolation and moving at a typical velocity of Vstar 100 km/s 
with respect to the external medium, the floor of interstellar den¬ 
sity value of ~ 1 cm“® will be limited by the ram pressure of the 
wind against the external medium, and will be set at a distance of 
( I Fryer, Rockefeller & Young||2006l) 


fwind = l-SpC 




lO-sjWeyr-i 

-0.5 


'^wind 


1 cm“3 


'^sta 


1000 kms“i 
1 


100 kms“i 


(13) 


Since the SN blast wave decelerates at large distances > rwind, it is 
safe to assume that it already probes the ISM medium, and not the 
wind from its progenitor. 

If the star is not in isolation, but part of a young stellar clus¬ 
ter, then taking the Galactic center clusters as an example, which 
contain 10-15 per cent of all Wolf-Rayet stars in the Galaxy, we 
can find a typical de nsity as the winds of severa l stars interact with 
each other (see, e.g., |Mimica & Gianniosl20llh . Let us take a clus¬ 
ter of 7?c ~ 1 pc size with N 100 0-stars, with t ypical wind s 
of Mwind 10 “®Mq yr“^ and Wwind 1,000 km/s dFigej2003) . 
The stellar separation is d ~ RcN~^^^. At the place were winds 
of different stars interact with each other, the density is typically 
of order 1 cm~®. This can be obtained by calculating the typical 
wind density at a distance ~ d, which is Mwind/(TTrd^iiwind) or, 
alternatively, by adding the winds of all stars in the cluster, which 
yields a total wind density of A^Mwind/(4:7ri?^r;wind), both expres¬ 
sions yield 1 cm~® density. Although the density as winds in¬ 
teract will not be constant, n ~ 1 cm“^ is a fair guess for the 
characteristic density of the ambient gafl In the unlikely case that 
this constant density floor is not reached and the wind-like medium 


^ The ISM may be clumpy (as revealed in, e.g.. Galactic remnants), 
and that makes the bl ast appearance much more rich (see, e.g., 
lObergaulinger et ^I20l4) . Since we are interested in the average proper- 


extends to hundreds of pc scales, then the SN blast wave will de¬ 
celerate > 5000 yr after the explosion for the typical parameters 
considered above. In this case, the strong and distinct radio SNR 
rebrightening discussed in this paper will not occur. 

In this paper we have taken a simplistic approach and com¬ 
pletely separated both components: the GRB and the SN ejecta. 
It is possible that both components are “connected” in the energy- 
velocity space, and there is a continuum of components between the 
relativistic GRB one and the non-relativistic SN one, which could 
be tested with future observations. However, this seems not to be 
the case for GRB 030329. For this burst, it appears that a single 
component is enough to model the GRB radio afterglow, and the 
presence of another component could have already potentially been 
observed. If a continuum of components exists between the GRB 
and the SN one, then energy would continuously be injected to the 
blast wave while it decelerates. This may also produce a distinct re¬ 
brightening in the light curve. For a constant density medium, as the 
energy injection transitions to the non-relativistic phase, the decay¬ 
ing radio light curve is expected to start rising as long as the energy 
is injected steeper than oc where /3 is the blast wave ve- 

locit}!^ [this statement is s omewhat dependent on the precise radio 
spectrum, see eq. (23) in [Barniol Duran et al.ll201^ . In addition, 
the rise in the light curve would also appear if the energy injec¬ 
tion is accompanied by a transition of the external medium from a 
wind-like to a constant density one. 

Nonetheless, a connection in energy-velocity space between 
the SN and GRB components is expected for low-luminosity 
GRBs (//GRB) in our sample, where the GRB em ission is likely 
produ c ed in the shock breakou t scenario (e.g. Kulkamiet_^ 
19981: iMatzner & McKeel 1 19991: iTan. Matzner & McKed l200ll: 


Campana et 20061 and recently. iNakar & Sarill2012h . Here, the 


explosion energy is deposited at the center of the progenitor and 
this drives a shock that crosses the star and accelerates at the stel¬ 
lar edge. This acceleration dicta tes a specific relation betw een the 
fast and slow moving material jMatzner & McKe^ll999t) . which 
is seen in regular SNe, but not in //GRB, where there seems to 
be more energy in the fast moving ejecta (e.g. , Soderberg et al.l 
2006a|^ Marguttiet^ 20131 : [Barniol Duran et al.ll2015l) . Recently. 
Barniol Duran et alj(l2015h modeled the < 1 yr radio afterglow data 
of //GRBs within the relativistic shock breakout model. While it is 
tempting to use their model (and their results) to predict the SNR 
radio emission, there are a few reasons why we have decided not 
to use it. First, as mentioned before, the < 1 yr afterglow data 
probes the external density at distances much closer to the source 
than the larger distances probed by the SN component. Secondly, 
since at face value, an extrapolation in energy-velocity space for 
these sources does not work, extrapolating this model all the way 
to /IsN and Ssn would be dangerous. Finally, iNakaJ ( 12015l) has 
recently suggested that the SN component is deposited at the cen¬ 
ter of the explosion, whereas the fast moving material is deposited 
by a (failed) GRB jet, thus decoupling both components. For these 
reasons, even for the //GRBs in our sample, we use the conserva¬ 
tive approach to clearly separate both components. Nevertheless, 
even in models where the interaction of the SN ejecta with the sur¬ 
rounding medium is invoked to explain the < 1 yr GRB radio data 


ties of a structure that is hardly resolved, a constant density medium may he 
a good first approximation. 

^ The energy-velocity of the GRB and SN co mponents in “regular” GRBs 
appears to he flatter than this condition (e.g., iMargutti et alj|2013h : how¬ 
ever, already GRB 030329 points out that the two components might not he 
connected, hut otherwise be well-separated. 
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fe.g.. lLi & Chevali^ll999l . for GRB 980425), a radio rebrighten¬ 
ing at late times will still occur. Here, the SN blast wave velocity 
is initially found to be > 0.5c, and the energy ~ 10®° erg, whereas 
> 100 times more energy is carried by the ~ O.lc-moving mate¬ 
rial, which will decelerate at much later times and can power a dis¬ 
tinct peak in the emission surpassing (long after the explosion) the 
emission of the GRB afterglow or other mildly relativistic compo¬ 
nents (unless the wind-like medium found for GRB 980425 extends 
to very large distances - see discussion above). 

For simplicity, we have not included the contribution from 
the GRB “counter jet” in our calculations: the second GRB jet 
that points away from us. Although there could be some confu¬ 
sion between the possible re-brightening of the counter jet and 
the SNR radio emission, the energy in the SN component is much 
larger than that of the counter jet, making the SNR emission much 
brighter. Also, the SNR radio emission is expected to peak at a 
few tens of years after the explosion, whereas the counter jet con¬ 
tribution t o the GRB radio after glow sh ould peak in a few years 
time-scale jGranot & Loeb||2003 . see also l Sironi & Gianniosl2013l 
and simulations bv ivan Eerten et alj|201(]|) . Even in the case of a 
wind medium, the GR B counter jet should peak before ~ 10 yrs 
jPe Colle et al.ll2oT^ . For the case of GRB 030329, there is no 
clear evidence yet for a re-brightening due t o the counter jet (e.g., 
Ivan der Horst et n]|2Q08l : lMesler et alj2012t) . 

There are several regular type Ibc SNe (not associated with 
GRBs) with upp er limits on their radio emission at late times af¬ 
ter the explosion iSoderberg et al.ll2006bl) . These SNe usually have 
kinetic energy of ~ 10®^ erg with velocity of 10,000 km/s. As 
a result of the low expansion velocity (of this particular compo¬ 
nent - see below), these SNe are expected to be fainter radio emit¬ 
ters. For these SNe, we c an also predict the expe cted SNR radio 
emission. In the sample of ISoderberg et alj l2Q06bll . SN1985F pro¬ 
vides the strongest constraint since it is a nearby source (7.7 Mpc) 
with a strong upper limit (37/iJy at 8.46 GHz, ~ 18 yrs after the 
explosion). This constrains the external density for this source to 
be <0.7 cm~®. Weaker density constraints can be determined for 
other SNe in this sample. On the other hand, typical Type Ibc SNe 
do show bright radio emission that peaks on ~ tens of days time- 
scale; however, this emission is produced by a O.lc-component 
with kinetic energy of 10“^® — 10“^® erg. 

In addition to predicting the SNR radio emission for GRBs 
with SNe association, we have also predicted it for GRBs with no 
optical SNe identification to very strong limits. The detection of a 
SNR radio emission would settle if a 1998bw-like event was truly 
present in these explosions. Similar predictions could be done for 
short GRBs at low redshifts, for which an accompanying SN is not 
expected and is not observed either to strong limits. For example, if 
a 1998bw-like event accompanied the short GRB 080905A, which 
is at a redshift similar to GRB 060614, then a similar SNR radio 
emission than the one calculated for GRB 060614 (see Fig. a just 
shifted in time by ~ 2 yr) would be present. Similar pre dictions 
could be made for other nearby short GRBs (see fig. 2 in lBerge^ 
l2014h . 


The detection of the SNR radio emission would help to con¬ 
strain the particle acceleration and magnetic field generation mech¬ 
anisms for < 0.1c shocks. It will also allow us to constrain the ex¬ 
ternal medium few pc scale from the center of the GRB explosion. 
Even if it turns out that the SNR radio emission is not observed, it 
will inform us that, for some reason, particles are not accelerated in 
< 0.1c shocks, and this would help guide current studies of particle 
acceleration models. We strongly encourage regular radio follow¬ 
ups of the locations of GRB 980425 and 030329. Any burst with 


2 < 0.2 that is more than one decade old also makes a promising 
target for radio follow-up. 
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